o 

^ : THE LARGE-SCALE CLUSTERING OF MASSIVE DARK MATTER HALOES 



^ 






c^ 



O 
O 



V. DESJACQUES 

Institute for Theoretical Physics, University of Zilrich 
Winterthurerstrasse 190, 8057 Zilrich, Switzerland 



\^ • The statistics of peaks of the initial, Gaussian density field can be used to interpret the 

r) I abundance and clustering of massive dark matter haloes. I discuss some recent theoretical 

results related to their clustering and its redshift evolution. Predictions from the peak model 

are qualitatively consistent with measurements of the linear bias of high mass haloes, which 

also show some evidence for a dependence on the halo mass M at fixed peak height v. The peak 

Q . approach also predicts a distinctive scale-dependence in the bias of haloes across the baryon 

)—{ ' acoustic feature, a measurement of which would provide strong support for its validity. For 

C/5 , 2a density peaks collapsing at z — 0.3, this residual scale-dependent bias is at the 5-10% level 

and should thus be within reach of very large simulations of structure formation. 



^ ' 1 Peaks in Gaussian random field 

in 

o 



The peak model introduced by^ assumes that dark matter haloes are associated with peaks 
of the initial (Gaussian) density field. Although dark matter haloes are the local maxima of 
lO ' the evolved mass distribution, there is a clear correspondence with initial density maxima for 

massive objects only. In the following, I will focus on the large-scale clustering properties of 
initial density peaks and show there is nontrivial scale-dependence both in the linear spatial and 
velocity bias. I will discuss some implications of these results. 
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H 2 First order biasing of initial density peaks 

FollowingLUj one usually smoothes the initial density fluctuations at redshift Zj ^ 1 with a filter 
of characteristic mass scale M before identifying local maxima of height i'. Even though density 
peaks form a well-behaved point process, the large-scale asymptotics of the 2-poin t co rrelation 
and pairwise velocity can be though of as arising from the continuous bias relationl^^ 

6np],{x) = bjM{x) -b(;A6Ai{x), (1) 

2 

Vpk(a;) = VMix) ^V6m{x) , (2) 

where (Jnpk and i^pk are the peak count-in-cell density and velocity, 5m and vm are the initial 
mass density and velocity field smoothed on scale M, and the (Lagrangian) bias parameters b^, 
and b(^ are 

bu{i^,li) = — i- V)' ^c(^'7i) = — (-j V) • (3) 

ctq V 1 - 7i / o"2 V 1 - 7i / 



Here, u = u{v) denotes the mean curvature of peaks of height v, 7i(M) = crf/(Too-2 and o"o, cJi 
and (72 are spectral moments which depend upon the shape of the hnear mass power spectrum. 
Note that h^ is strictly positive, whereas hy can be positive or negative. In Fourier space, 
wavemodes of the peak number density (5npk(fc) can be obtained by multiplying (5m (^) with 
(here and henceforth, I will omit the dependence on u and 71 for brevity) 

hp^{k) = K + h(^k'' . (4) 

This defines the spatial peak bias at the first order. In practice, the peak-background split 

o 
approach, which is based on count-in-cells statistics, can also be used to estimate by. In this 

regards, the linear Lagrangian bias h^, predicted by the peak model is exactly the same as that 
returned by the peak-background split argument*^. 

The peak velocity Vp\^{x) as defined in Eq.([2]) is consistent with the assumption that initial 
density peaks move locally with the dark matter. However, the 3-dimensional velocity dispersion 
of peaks is smaller than that of the mass (t_i, i.e. o"^ j, = a'^_i (1 — 7g) with 70 = o"o/o"_iO"i, 
because large-scale flows are more likely to be directed towards peaks than to be oriented ran- 
domly . As shown inl^, this leads to a fc-dependence of the peak velocities as can been seen 
upon taking the divergence of Vp]^{x) and Fourier transforming it, 

0pk(fc) = ( 1 - ^ kA W{k, M) e{k) = K,,{k) eM{k) , (5) 

where 9 = V ■ v \s the mass velocity divergence and W{k, M) is the Fourier transform of the 
filter. This defines the statistical velocity bias 6vei(fc)- Note that h^eiik) does not depend on u 
and, for the highest peaks, remains scale-dependent even though the spatial bias hp^{k) has no 
^-dependence in this limit. 

3 Redshift evolution of the peak correlation 

Pairwise motions induced by gravitational instabilities will distort the primeval peak correlation. 
The gravitational evolution of the correlation of initial density peaks can be addressed with the 
Zel'dovich ansatz^L^, assuming they behave like test particles moving with the dark matter. 
In this first order approximation, the gravitationally-evolved peak correlation ^pk(^, -s) is the 
Fourier transform of the peak power spectrum'^ 

Ppk(A:, z) = G\k, z) [K,,{k) + 6pk(A;, z)f Puih 0) , (6) 

where hp]^{k,z) = D{zi) / D{z)hp]^{k) and the function 

is a damping term induced by velocity diffusion. It is similar to the propagator Gs{k,z) intro- 
duced in"^, although the latter involves the matter velocity dispersion o"_i. The first term in the 
square bracket reflects the fact that peaks stream towards (or move apart from) each other in 
high (low) density environments, but this effect is A;-dependent owing to the statistical velocity 
bias. Therefore, the Eulerian and Lagrangian linear bias parameters are related according to 



D(z) ^ D(z] ai 



h''A^)^l + ^^hMl 6^(z) ^ i^6c(z.) - 3 • (8) 



The first relation is the usual formula for the Eulerian, linear scale-independent bias^^. The 
second relation shows that 6? approaches —CTQ/al with time. 
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Figure 1: Large-scale bias of dark matter haloes identified with a FOF finder of linking length b = 0.2 (left 

panel) and a SO finder with redshift-dependent overdensity threshold (right panel). Circles and triangles refer to 

halo samples whose average mass is ~ 1.3 and 5 x 10^'^ Mq/Zi, respectively. The dotted and dashed curves are 

formulae based on the excursion set theory, whereas the solid curves are fits motivated by the peak model. 



4 The large-scale bias of dark matter haloes 



The large-scale bias contains important information on the abundance and clustering of biased 
tracers of the density field. To compare theoretical expectations with measurements of dark 
matter halo bias, I will assume that peaks of height i' = 6c/cro{R,z) identified in the initial, 
smoothed density field 5m are associated with objects of mass M collapsing at redshift z. 

The peak model predicts that, for moderate peak height, b^ is significantly smaller than the 
value 1 + u'^/Sc derived for thresholded regions^^ due to the correlation between the peak height 
and the peak curvature^. However, in the limit u ^ 1, b^{i') ~ 1 + (i/^ — 3)/6c which shows 
that the evolved linear bias of initial density peaks of height v indeed converges towards the 
prediction or^'^. This should be compared to well-known expressions derived from the extended 
Press-Schechter formalism which, in the same limit, evaluate to &mw(^) = 1 + (^^ ~ ^)/5c^^ 
and b^r^{u) ~ 1 + {av'^ — l)/6c'^^- In the latter case, a = 0.75 follows from normalising the 
Sheth-Tormen mass function to N-body simulations. Note that, whereas fe^w ^^'^ ^ST depend 
only upon the peak height, b^ is a function of both v and M (through 71 (M)). 

In Fig. [H these various predictions are compared with measurements of the linear bias of 
massive haloes extracted from numerical simulations of structure formation -. Error bars show 
the scatter among various realisations. The measured halo bias appears to depart fron i the 
Sheth-Tormen scaling at large v, in agreement with recent measurements of the halo bias ^^ i 
Furthermore, the data shows evidence for a dependence on M, but the exact magnitude of the 
effect is sensitive to the halo finder. Because the best choice of filter is a matter of debate, I 
treat 71 as a free parameter and show 6^(z^, 71) for 71 = 0.4 and 0.5 (a Gaussian filter yields 
71 ~ 0.65 for the mass range considered), which provide a reasonably good fit to the bias of 
^2cj haloes. Note that the peak expression b^ is also found to match the bias of massive haloes 

in scale-free cosmologies rather welll^. 



5 Peak biasing and the baryon acoustic oscillation 

Having checked that the peak model predicts a large-scale halo bias b^{z) consistent with simu- 
lations, I consider now the impact of the scale-dependent piece bf{z)k'^. The presence of such a 
term amplifies the contrast of the baryon acoustic oscillation (BAO) in the correlation of initial 
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Figure 2: Left: Redshift evolution of tiie baryon acoustic oscillation in the correlation of initial, 2a density peaks 

as predicted by Eq.®. Results are shown at redshift z = 100, 5, 2, 1 0.5 and 0.3 (curves from bottom to top). 

Middle: Same as left panel but for "linear tracers", for which the correlation simply is b^{z)^ times the evolved 

mass correlation Right: The ratio diverges at r ~ 130 h~^Mpc because zero-crossings do not coincide 

density peaks relative to that in the hnear theory correlation l^. Eq.([6|) can be used to estimate 
how much of this effect survives at virialization redshift (a more realistic calculation should 
include the mode-coupling power). 

To emphasise the effect of bf{z)k'^, Fig. [2] compares the redshift evolution of the large-scale, 
2-point correlation ^p^ of initial density peaks (left) with that of "linear tracers", .^n, for which 
Pit{k,z) = G'^{k,z)[b^{z)]'^ PM{k,0) (middle). The right panel displays the ratio between the 
two correlations. Results are shown for 2a density peaks collapsing at Zc = 0.3 and identified 
on a mass scale 5 x 10^^ Mq/Zi with a Gaussian filter. The relative amplification of the BAO 
contrast in Cpk(^, -^i) induces a scale-dependence in the bias that decays with time owing to the 
smearing from velocity dispersion. At the collapse redshift however, the model predicts residual 
scale-dependence across the BAO feature at the 5-10% level (right), a measurement of which in 
numerical simulations would provide strong support for the validity of the peak approach. 
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